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Abstract 

We consider the pattern formation problem in coupled identical systems 
after the global synchronized state becomes unstable. Based on analytical 
results relating the coupling strengths and the instability of each spatial mode 
(pattern) we show that these spatial patterns can be selectively realized by 
varying the coupling strengths along different paths in the parameter space. 
Furthermore, we discuss the important role of the synchronized state (fixed 
point versus chaotic attractor) in modulating the temporal dynamics of the 
spatial patterns. 
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Introduction. In a reaction-diffusion system, at the Turing bifurcation point the global 
equilibrium state becomes unstable and an inhomogeneous spatial pattern known as a Turing 
pattern [1] emerges. Following Turing's classic work pattern formation in reaction-diffusion 
systems has become a major topic of investigation both theoretically and experimentally 
[2,3]. 

The idea in Turing's original formulation is not limited to just reaction-diffusion systems 
with a single diffusion constant. Networks where the diffusion is characterized by a more 
complex function and networks where the coupling is not diffusive can be studied in a 
similar fashion. This makes the idea of pattern formation applicable to diverse areas of 
science and engineering [4-11]. Furthermore, the globally synchronized state prior to the 
formation of spatial patterns need not be a fixed point as in Turing's original case. One 
can consider synchronized limit cycles and chaotic attractors as the bifurcating state. When 
the synchronized state becomes unstable under parameter variation, the dynamics may exit 
the synchronization manifold along certain eigen-mode, giving rise to a spatial pattern. We 
refer to these patterns as Generalized Turing Patterns (GTPs). 

GTPs have been the subject of several recent publications [10]. These papers mainly rely 
on numerical techniques to obtain the threshold of instabilities. In this letter we focus on 
the analytical treatment of the pattern formation and pattern selection. To do so we restrict 
ourselves to networks in which each node is coupled to its P nearest neighbors. Under the 
assumption of periodic boundary condition and translational invariance the eigen-modes are 
the spatial Fourier modes. Our main results are as follows. First, we derive explicit analytical 
expressions defining the stability region in the parameter space spanned by the coupling 
strengths. Second, we demonstrate that, by following different paths in the parameter 
space, different spatiotemporal patterns can be selectively realized. Third, although the 
spatial patterns are the same for a given coupling matrix, we show that, depending on 
whether we have a synchronized fixed point or chaotic attractor, the temporal evolution of 
the patterns is either constant in time or modulated in an on-off intermittent fashion. We 
study networks on both one dimensional (1-d) and two dimensional (2-d) lattices. Only 
coupled maps are considered here. Similar results apply to coupled differential equations. 

One dimensional lattices. Consider the following 1-d coupled map lattice where each 
node is coupled to its P nearest neighbors in a general diffusive manner: 

x,(n + 1) = f (x») + -^£> p [f (Wn)) + f (x,_») - 2f (x»)] , j = 1, 2, . . . , L. 

(1) 

Here the M-dimensional vector x.,- describes the map at the jth node of the lattice, a p 's 
denote the coupling strengths and n is the time index. We impose the periodic boundary 
condition: x J+ L(n) = Xj(n). Since the coupling strengths don't depend on the spatial 
position index j we have translational invariance. The individual (uncoupled) map is Xj (n + 
l)=f(x»). 

Clearly the synchronized state x.,- = x, V j is a solution of the above system. To study 
its stability we linearize Eq. (1) around the synchronized state x to obtain: 

1 p 

Zj(n + 1) = J(x(n))z j (n) + J(x(n)) — ^ a p [z J+p (n) + z j _ p (n) - 2z j (n)}, (2) 

Ar p=i 



where Zj denotes the jth map's deviations from x and J is the M x M Jacobian matrix 
corresponding to f. Consider the discrete Fourier transform of zf 



1 L 



= tE exp(-i27rj7/L)zj(n). (3) 



Performing the discrete Fourier transform on the linearized equation we get 



1 p 



»7,(ra + 1) = J(x(ra)) 
Simplifying we finally have 

rj^n + 1) = J(x(ra)) 



*7i( n ) + 7T5 J2 a v (exp(i2Trpl/L) + exp(-i2irpl / L) - 2) 17,(71) 
A ^ P =i 



• (4) 



2 p 

1 - — Ea p sin 2 (7rp//L) 



m{n). (5) 



We start by examining the case where the individual maps are chaotic and x corresponds 
to the synchronized chaotic state. The synchronized state is stable if all the transverse 
Lyapunov exponents are negative [11,12]. Computing the Lyapunov exponents /!« of the 
above linearized equation we get 

2 p 

= /i l + ln|l--Ea p sin 2 (7rp//L)|, i = l,2,...,M, (6) 
F P =i 

where h^s are the Lyapunov exponents for the individual map ordered as h± > h 2 > • • • > 
h M - Note that the / = Fourier mode corresponds to the synchronized chaotic state where 
the Lyapunov exponents of the coupled system are given by the Lyapunov exponents of 
the individual system. For each I 7^ value, gives the largest transverse Lyapunov 

exponent. Therefore, stability of the synchronized chaotic state is ensured if < for 
all I 7^ 0. By symmetry of the Fourier modes, we only need to consider < for 

I — 1, 2, . . . L/2 (if L is odd, we take (L — l)/2). Thus we have the stability conditions 

2 p 

|l-pE« P sin 2 (7rp//L)| < exp(-/i 1 ), 1 — 1,2,... L/2 or (L — l)/2. (7) 
F P =i 

This set of inequalities defines a stability region in the parameter space spanned by the 
Op's. By selecting a given Fourier mode and choosing a suitable path in the parameter space 
we can realize the corresponding GTP. Note that, if one considers only the nearest neighbor 
diffusive coupling (P = 1), the parameter space is one dimensional and at most two GTPs 
can be excited by varying the coupling strength even though there are many possible modes. 
By enlarging the parameter space we obtain much greater variety in terms of GTPs that can 
be realized. With recent progress [13] in obtaining stability conditions for general coupling 
matrices, we are currently investigating whether even more general Turing patterns can be 
obtained. 

As a numerical example we consider coupled logistic maps in the chaotic regime where 
f(x) — 1 — ax 2 with a = 1.9. The maximum Lyapunov exponent hi is 0.549. For simplicity, 



we restrict ourselves to L = 5 and P = 2. The stability conditions for the synchronized 
chaotic state are: 

-0.578 < 1 -ax sin 2 (7r//5) -a 2 sin 2 (27r//5) < 0.578, Z = l,2 (8) 

which give the stability region marked black in the parameter plane [Fig. 1(a)]. We call the 
I = 1 mode the long wavelength (LW) pattern and the / = 2 mode the short wavelength 
(SW) pattern. The arrows indicate the paths along which either of the two patterns can be 
selected. 

The main frame in Fig. 1(b) shows the temporal dynamics of the long wavelength 
pattern for a x = 0.96 and a 2 = 0.1. Here deviations from the synchronization manifold is 
approximated by 

L 

Zj[n) =x j (n)-^2x j (n)/L 
j'=i 

with L — 5. To facilitate visualization, at each time step n, a continuous function is splined 
through the six discrete nodes: Zi(n), ^(n), • • •, 25 (n), and Zq(ti) = Z\{n). Furthermore, to 
overcome the distortion due to the two opposite phases of a pattern, we monitor the deviation 
at a given node and multiple the deviations at every node by a minus sign whenever the 
deviation at the monitored node becomes negative. 

Since the bifurcation undergone by the system at the boundary of the stability region 
is the blow-out bifurcation and there is only one attractor prior to the bifurcation, the 
dynamics in this case is referred to as on-off intermittency [12,14]. The temporal evolution 
of the deviations at a typical node is given by the curve to the left of the main pattern frame. 
Its bursting behavior is characteristic of on-off intermittency. The GTP itself is given at the 
bottom of Fig. 1(b). 

For ai = 0.04 and a 2 = 1.1 we observe the short wavelength pattern in Fig. 1(c). The 
same visualization methods are used to make this figure. 

To understand how the synchronized state shapes the temporal properties of the GTP 
after desynchronization we consider a case where the synchronized state is a fixed point. 
Again we use the coupled logistic maps and choose a = 0.5. The fixed point is x = 0.73 and 
the lyapunov exponent is hi = —0.31. Still letting L = 5, P = 2, we get the stability region 
defined by 

-1.36 < 1 - aisin 2 (7r//5) — a 2 sin (2-7T//5) < 1.36, 1 — 1,2 (9) 

shown in Fig. 2(a). Following the arrows we can realize either the long or the short wave- 
length patterns. Figure 2(b) gives the long wavelength pattern for a± = 0.5 and a 2 = 2.5 and 
Fig. 2(c) gives the short wavelength pattern for ai = 2.5 and a 2 = 0.5. The same methods 
of plotting as that used for Figure 1 are used here. Comparing Figs. 2(b) and (c) with Figs. 
1(b) and (c) we see the same spatial patterns but different temporal behaviors. 

For the example in Fig. 2 the final GTPs are the new fixed points displayed as a function 
of the space coordinate. Predicting the exact location of these new fixed points requires the 
nonlinear terms dropped in the linear stability analysis. Although it is often the case that 
the spatial functions underlying the new fixed points agree with the respective linear eigen- 
modes this is by no means a guaranteed fact [2] . On the other hand, when the synchronized 



state is chaotic, linear analysis will govern the temporal evolution whenever the phase space 
trajectory comes back to near the synchronization manifold. 

Two dimensional lattices. Finally, we consider a 2-d coupled map lattice given by 

1 p 

Xj )fc (n + 1) = f(xj- fc (n)) + — ]T {a p [f(x j+M (n)) + f(xj_ Pifc (ra))] 

Zjr p=i 

+6 P [f (x j)fc+p (n)) + f (x jifc _ p (n))] - 2(a p + 6 p )f (x i)fc (n))} , j, k = 1, 2, . . . , L. (10) 
Linearizing around the synchronized state x J)fc = x, V j, fc, we get 

1 p 

Zj,k( n + !) = J( x W) z j,fc(«) + J ( X W)^E W z i+P,fc( n ) + z i-P,fcH] + 

Zjr p=i 

bp [ z i,fe+p (n) + Zj jfc _ p (n)] - 2(a p + 6 p )z iifc (n)} , (11) 

where z denotes the deviation from the synchronized manifold. Applying the 2-d discrete 
Fourier transformation given by 

y L L 

Vi,m( n ) = y^^^exp(-i2vrj7/L)exp(-i27rA;m/L)z iifc (n), (12) 

L j=l k=l 

to the linearized equation and simplifying we finally obtain 



m,m(n+l) = J(x(n)) 



2 p 

1 — — ^2 (a p sin 2 (7r pi / L) + b p sin 2 (irpm / L)^J 
P P =i 



Vi, m (n). (13) 



Proceeding as before, the stability conditions for synchronized state are given by 
2 p 

|1 — — (a p sin 2 (tt pi /L) + b p sin 2 (npm/L)^ \ < exp(—hi), 
P P =i 

l,m = 0,1,2,... L/2 or (L - l)/2, (14) 

with l — m — excluded (as it corresponds to the synchronized manifold). 

As in the 1-d lattice case, GTPs emerge when the synchronized state loses its stability. 
This happens when the parameters are varied across the stability boundary. We can select a 
particular (I, m) mode and have it realized by appropriately choosing the coupling strengths 
a p and b p (p = 1, 2, . . . , P). We will again illustrate this with logistic maps coupled together 
on a 2-d lattice. We set a = 1.5. Then hi = 0.231. Choose L = 5 and P = 1 for simplicity. 
The stability conditions are 

-0.79 < 1 - —sin 2 (7r//5) - —sin 2 (7rm/5) < 0.79, /,m = 0,l,2, (15) 

with I — m — excluded. Figure 3(a) depicts the stability region (black). Different edges of 
the region are the instability thresholds of different modes. It is obvious that only the (0,1), 
(1,0) and (2,2) modes can be observed given two tuning parameters. By including additional 
coupling strengths not considered here one can potentially make the (1,1), (1,2) and (2,1) 
modes observable. Figure 3(b) shows a snapshot of deviations from the synchronization 



manifold at different nodes for a\ = 0.5 and b\ = 0.5. As before a continuous function is 
fitted through all the nodes. Although not perfect the essential features of a (2,2) pattern is 
quite clearly visible in the snapshot. The other two patterns (0,1) and (1,0) can be similarly 
realized by choosing proper parameters with respect to the boundary of the stability region 
in Fig. 3(a). 
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FIGURE CAPTIONS 



Figure 1: Pattern selection from the synchronized chaotic state in a 1-d map lattice 
(P = 2). In (a), the region of stable synchronization (black area) and distinct pattern 
selection directions are shown. In (b), temporal evolution of the long wavelength 
pattern is given with a\ = 0.96, 02 = 0.1. In (c), temporal evolution of the short 
wavelength pattern with a± = 0.04, a 2 = 1.1 is given. 

Figure 2: Pattern selection from the synchronized equilibrium state in a 1-d map lattice 
(P = 2). In (a), the region of stable synchronization (black area) and distinct pattern 
selection directions are shown. In (b), the long wavelength pattern evolving as a fixed 
point in time is given with ai = 0.5, a 2 = 2.2. In (c), the short wavelength pattern 
with ai = 2.2, a 2 = 0.5 is shown. 

Figure 3: Pattern selection in a 2-d map lattice (P = 1). In (a), the region of sta- 
ble synchronization (black area) and the directions of selecting some patterns in the 
parameter space are given. In (b), a (2,2) pattern is shown with ai = 0.5, 61 = 0.5. 
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